Recycled Noise Rectification: A Dumb Maxwell's Daemon 
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The one dimensional motion of a massless Brownian particle on a symmetric periodic substrate 
can be rectified by re-injecting its driving noise through a realistic recycling procedure. If the 
recycled noise is multiplicatively coupled to the substrate, the ensuing feed-back system works like 
a passive Maxwell's daemon, capable of inducing a net current that depends on both the delay 
and the autocorrelation times of the noise signals. Extensive numerical simulations show that the 
underlying rectification mechanism is a resonant nonlinear effect: The observed currents can be 
optimized for an appropriate choice of the recycling parameters with immediate application to the 
design of nanodevices for particle transport. 

PACS numbers: 05.60.-k, 05.40.-a, 74.50. +r 



I. INTRODUCTION 

When control or spurious signals, cither periodic or 
noisy, are injected into an extended system, they can cou- 
ple a variable of interest x(t) through different channels, 
additively or multiplicatively, alike. While being trans- 
mitted across the system components, an input signal 
may undergo 0, Q : 

(a) time delay, due to the combination of diverse prop- 
agation or transduction mechanisms. As a result, an in- 
put signal rj(t) can split into two or more forcing signals 

£i(t) acting on x(t), each with a time delay i.e. 
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Let us consider, for instance, an electric signal acting 
on a particle of coordinate x(t) in a narrow channel after 
propagating e.g. through two different nondispersive me- 
dia, as illustrated in Fig. Qfa). Regardless of the nature 
of the input signal, a wave train or random pulses, the 
time delay Td of Ba(t) relative to E\(t) depends of the 
electromechanical properties of the medium surround- 
ing the channel; moreover, the two drives are likely to 
couple the particle differently: E\{t) pulls the particle 
horizontally and, therefore, additively, while E 2 (t) prop- 
agates normally to the channel, thus modulating multi- 
plicatively its effective substrate potential; 

(b) time correlation, due to the finite response time 
of the system components. The relevant response, or 
transfer functions modulate the amplitude of the signals 
effectively coupled to the variable of interest x(t) 
and, more importantly, determine a finite autocorrela- 
tion/decay time Tr^ for each depending on its cou- 
pling channel [l|,|s3- 

Figure ^b) shows the ideal case of a particle diffusing 
through a one dimensional (ID) pore, or channel in a ho- 
mogeneous dielectric medium with relaxation time r c ; the 
sample is placed between the plates of a capacitor sub- 
jected to a variable voltage. The particle is thus directly 
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FIG. 1: (Color online) Effective action of an external field 
of force on a Brownian particle in a narrow channel: (a) The 
external field splits into a direct, E\(t), and a transmitted 
component, E<z{t), that couple to the particle with relative 
delay t c ; (b) The field couples to the particle both directly, 
Ei(t), and modifying its environment through the local field 
£?2 (t) . In both cases E2(t) modulates the channel geometry 
(multiplicative coupling) . 



coupled to the time-dependent electric field Ei(t); how- 
ever, its mobility in the x direction changes in time as the 
channel cross-section is modulated by the electrostriction 
effects induced by the capacitor. The potential energy of 
the particle along its path is thus a function of the local 
field E^(t), which, in turn, can be regarded as a dielectric 
response to the input E\ (t) . As long as the frequency de- 
pendence of t c can be neglected (low frequency regime), 
amplitude and phase-shift of each spectral component of 
Ei (t) can be explictly computed in linear response theory 
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1] . If Ei (f ) is a random signal with correlation time r c , 
then E 2 (f ) is an autocorrelated noise with time constant 

ri 2 -* = t c + Tc 3], while the primary, Ei(t), and the 
secondary signal, E-2.it), are crosscorrelated with effec- 
tive time constant r c . If E±(t) is an ac field with angular 
frequency w, the signal E2{t) modulating the channel de- 
velops a time delay r c ~ arctan(a;T c )/a> [j]. 

In this paper we work with stationary Gaussian noise 
sources and restrict ourselves to the case of two distinct 
coupling channels, only, so that a given noise source r)(t) 
splits into two colored driving noises, 77(f) — > £2(^)1 
with correlation times ri 1 ^ and Tc , and delay times 
and . To simplify our notation we further assume 
that: (i) r c (1) = r c (2) = r c > 0; (ii) rf 3 - = r d > 0, 
which for a stationary 77(f) is equivalent to setting r^p = 

and r^ 2) = Td. 

For the purpose of numerical simulation, conditions (i) 
and (ii) are satisfied by noises of the form II) with 77(f) 
a Gaussian noise with autocorrelation function 3] 

(77(f)77(0)> = e-^/r c . (1.2) 

All correlation functions (£i(f)£j(0)), with i,j — 1, 2, can 
be readily expressed in terms of the autocorrelation func- 
tion (|I.2(I . According to this notation, for practical pur- 
poses we term £i(f) the primary or source noise, and 
^(f) the secondary or recycled noise [with no reference 
to '/(/ : • 

Physical systems driven by delayed correlated noises 
are rather common in nature, typical examples being the 
propagation of charge density waves HH, the migration 
of both pointlike and linear defects in crystalline materi- 
als , the transport of nano-particles in biological § and 
artificial channels [t| , the manipulation of vortex lines in 
superconducting devices 0] and colloidal particles along 
ID tracks to mention but a few. 

In most models discussed in the literature, the input 
signal 77(f) is time periodic and so are the two (or more) 
driving terms £j (f ) |q.p2.p3j| : A given phase-lag between 
two additive signals I6|.|l2| or between an additive and a 
multiplicative signal |l3lll4| may breach the spatial sym- 
metry of the underlying x dynamics, thus inducing a net 
current (x) . [Note that the average of (x) over a uniform 
distribution of the phase-lag vanishes (spontaneous sym- 
metry breaking mechanism).] Here, we consider similar 
models with the difference that the ac drives are replaced 
by the noises £j(f) introduced above. 

Our models with two (or more) delayed noises should 
not be mistaken for the nonlinear systems with delay re- 
cently investigated by several groups in the context of 
stochastic resonance 0|: These authors proposed 
to control the dynamics of x(f) by utilizing appropri- 
ate functions of x(t — T4) as feedback terms. This is 
a well-established control technique in physics and elec- 
trical engineering 0. Extended experimental appara- 
tuses, where both types of delays (from recycling and 
feedback) must be taken into account, are, for instance, 



the gravitational-wave interferometers, like the VIRGO 
detector pj|. Here, an external signal 77(f), e.g., a seis- 
mic disturbance, enters the antenna by creeping through 
its mechanical suspensions and dampers and eventually 
combines with the intrinsic electronic and photonic noises 
of the apparatus, so that the detection signal x(t), cor- 
responding to the mirror displacement induced by the 
gravitational signal, is additively and multiplicatively af- 
fected by 77(f) at different times. Moreover, the control 
loops that maintain the alignment and the locking of the 
interferometer, make use of sophisticated feedback tech- 
niques also involving x(t). 

In all examples mentioned so far, the multiple action of 
an external disturbance is often regarded as a nuisance; 
on modeling the system of interest one simply takes no- 
tice of its presence and tries to minimize its impact on 
the system response. In our approach we took a more 
"pro-active" stance, namely, we propose to tap the pri- 
mary noise source £i(f) and re-inject the recycled noise 
signal £ 2 (f) into the system so as to control its response. 
At variance with an earlier attempt |19| | , we consider here 
the motion of a massless Brownian particle diffusing on a 
symmetric ID periodic substrate; for appropriate choices 
of the time constants r c and t ( i the particle drifts sidewise 
with net velocity (x). Such a mechanism can be viewed 
as an automated Maxwell's daemon, where the rectifica- 
tion of the primary noise signal £1 is achieved through 
a pre-assigned recycling protocol: Lacking the dexterity 
originally assumed by Maxwell, our "dumb" stochastic 
daemon performs its chore " in average" and with reduced 
(but dependable) efficiency. Note that rectification of a 
primary noise by nonlinear interference with a recycled 
image of itself can hardly be assimilated to a ratchet 
mechanism |2jj . as here no spatial substrate asymmetry 
is required. 

This paper is organized as follows: In Sec. II the pri- 
mary, £1, and the secondary noise, £2? are coupled ad- 
ditively to an overdamped variable x bound to either a 
parabolic well (Ornstein-Uhlenbeck process) or a quartic 
double well (Kramers problem) p|. In Sec. Ill we simu- 
late numerically a massless Brownian particle diffusing, 
under the action of an additive primary noise £1, in a ID 
cosine potential modulated in amplitude by the recycled 
noisy term £2- The particle rectification current (x) at- 
tains an optimal intensity in an appropriate range of the 
time constants t c and r^. In Sec. IV we discuss the de- 
pendence of (x) on the remaining simulation parameters 
in view of the potential applications of our rectification 
scheme to the design of nano-devices for particle trans- 
port. 



II. NOISE RECYCLING 

The technique of noise recycling is well-known in laser 
interferometry |17|: One extracts a signal from the fluc- 
tuation source to be suppressed and re-injects it into the 
system after appropriate manipulation; the process in- 
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volves a certain number of steps, like analogue or digi- 
tal acquisition, usage of filters and actuators, etc, which 
eventually cause a delay of the feedback with re- 

spect to the primary signal £1 (t) . 

In order to familiarize the reader with the interplay 
of two delayed noises, we first consider the Ornstein- 
Uhlenbeck process 



■6(*)±&(t), 



(III) 



where £i(t) = VQir](t) and £ 2 (i) = VQzvit - T d), with 
T d > 0, and rj(t) denotes a zero-mean, white Gaussian 
noise with autocorrelation function {r)(i)rj(0)) — 28 (t). 
The stochastic differential equation (SDE) (|II.1|) can be 
treated analytically by mens of standard techniques 0,0- 
In particular, 



Q1 + Q2 



± 2 e 



a 



(1 - er 2at ) 



{11.2) 



e- 2a{t - Td) ]Q(t - r d ) 



where O(x) is the Heaviside step- function and asymptot- 
ically 



lim ([x(t) - x(0)] 2 ) 



(II.3) 



= -[{Q 1 + Q2)±2^Q 2 e- aTd ]. 



Analogously, one can compute the autocorrelation func- 
tion 



C{t) = lim (x(t + t)x(t)) 

T — >OC 



VQIQ 



2 e - at 



VQIQ~ 2 



± 



(V QT 



+ e 



(II.4) 

— a(|t— Td| — t 



with r d > 0. Note that C(-i) = C(i) and C(0) = (a; 2 ) > 
0, as it should. 

In Fig. [S]we compare numerical simulation results with 
our predictions 1|II.3J| . inset, and l|II.4|l . main panel: in 
all cases the corresponding data sets are indistinguish- 
able. Notice that sampling the ^-autocorrelation only 
for t > Td, thus overlooking the possibility that the pro- 
cess can be driven by two delayed noises, would lead to 
overestimating (x 2 ), 

, 1 



[(Qi + Q 2 ) ± cosh(ar d )], (II.5) 



with respect to the correct value (|II.3ll . Generalizations 
of the process Ijll.lf) to account for colored noises, with 
exponential autocorrelation functions of the type in Eq. 
(jl.2jl . and/or inertial effects, with x replaced by x + 75c, 
can also be treated analytically. Harmonic analysis Q 
yields straightforward, though cumbersome expressions 
for C(t). 

More suggestive is the nonlinear case of a massless par- 
ticle bound to a quadratic double well potential V(x) 
-ax 1 1 a + bx A /A 



x = ax 



•te 3 + 6(t)±6(*) 1 



(11.6) 




FIG. 2: (Color online) Autocorrelation function C(t) of the 
process HLH for Qi = Q2 = 2 and a = 2; the recycled noise 
£2 is added to (filled circles) or subtracted (empty circles) 
from the primary noise £1 after a delay time Td = 1- The 
solid curves represent the corresponding analytical predictions 
111.41 . Inset: C(0) = (a; 2 ) vs. r^; notation and the remaining 
simulation parameters are as in the main panel. 



with £1 and £2 defined as above, and a, b > 0. Sub- 
jected to the random kicks by the fluctuation sources, the 
overdamped Brownian particle hops between two bistable 
states centered at ± y/a/b and separated by the potential 
barrier AV = a 2 /46. According to Kramers' rate theory, 
the relevant hopping time 0, Q in the low noise regime 
is 



iry/2 



exp 



AV 



D±(r d ) 



(II.7) 



= (VQi + VQ2) 2 and 



where D±(rd) is the effective noise intensity of ± 
^2(^)1 t° be determined next. 

For zero delay, r d = 0, D_| 
D_ =0; for extended delays, r d — > 00, the two noises 
get completely uncorrelated, so that D± = Qi + 
As a consequence, the simulation results of Figs. E2a) 
and (b) reveal the following limiting behaviors of Tj^: (i) 
For r d — > co, both T^(r d ) tend to one limit, (|II.7|I . with 
D± = Qi +Q2 - perturbation corrections to Tj^(oo) may 
be computed in powers of D±/AV |]|; (ii) On subtracting 
the recycled from the primary noise, the Kramers' time 
Tft(T d ) tends to diverge; the effect gets amplified at small 
T d , when the destructive interference of the two noises 
is the most effective; (iii) When adding £1 and £2 with 
T d — > 0, the Kramers' formula l|II.7|) still applies but for 
D+ = {VQl + VQi? ■ Of course, T+(Q) < T|(oc). 

The T d dependence of T K (Td) for intermediate delays 
T d ~ ^k(°°) i s more interesting: The curves T^(T d ) 
exhibit a broad peak, whereas the curves T^(r d ) develop 
a shallow dip within approximately the same T d interval. 
This is a manifestation of the so-called resonant activa- 
tion |2(il | . Introducing a delay T d when superimposing two 
noises £1 and £ 2 with the same sign (+ sign), hinders the 
hopping dynamics for T d <C Tj^(oo); indeed, increasing 
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FIG. 3: (Color online) Activation process in the bistable 
potential V(x) = -ax 2 /2 + bx 4 /A driven by £i(£)±£ 2 (i): (a), 
(b) T£ (filled symbols) and (empty symbols) vs. Td for 
t c = and Qi = Qi\ (c) Probability density P(x) for r c = 0, 
Qi = Q-2 = 0.025, and = (circles), 0.3 (squares), and 100 
(triangles). 

Td decorrelates £i and £2 in time, i.e., diminishes the ef- 
fective intensity D + (rd)- This effect goes through a max- 
imum for Td ~ ^T^Td): A strong kick of £1 capable of 
making the particle reach the top of the barrier at x = 0, 
will be counter-balanced in average by an opposite kick of 
£2, which prevents the particle from falling into the other 
well. At even larger r^, the action of the two noises grows 
totally uncorrelated and T^(rd) approaches T^(oo) from 
above. At variance, combining £1 and £2 with opposite 
signs (— sign) produces a synchronization effect; for in- 
termediate delays Td hopping occurs after shorter waiting 
times, hence the dip in the T^(Td) curves. 

On lowering Qi and Q2 the Kramers' time increases 
exponentially; due to the prolonged sojourn at the well 
bottoms, memory effects in the particle dynamics are 
suppressed and so is the resonant behavior of T^(jd)- Fi- 
nally, we stress that the x probability densities P(x) obey 
the Boltzmann law P(x) oc e~ v ( x '' D± only for Td — ► 
and Td — > 00. For intermediate Td values the reduction of 
the dynamics (|II.5(I to a ID stationary process is ques- 
tionable 0, especially in the vicinity of the barrier; in 
other words, fitting the simulation P(x) curves of Fig. 



IHfc) by means of the Boltzmann law does not rigorously 
define D± in Eq. H11.7|) ~ but rather allows to extract a 
working estimate of it. 

Note that laser systems with o pto -electronic feedbacks, 
like those described in Refs. |16lll8| . are already capable 
of experimentally investigating the Kramers' rate mech- 
anism driven by delayed noises. 

III. A DUMB MAXWELL'S DAEMON 

At variance with the models of Sec. II, the recycled 
noise £ 2 can be re-injected so as to modulate the sub- 
strate on which the Brownian particle diffuses subjected 
to the primary noise £1 . In Ref . ^| we investigated the 
confined process 

x = a[l + &(t)]x - bx 3 + £i(t), (III.l) 

with £1 and £2 defined as in Sec. II and V(x) = —ax 2 /a+ 
bx 4 /4: with a,b > 0. The SDE i|III.l(l is not symmetric 
under reflection x — ► —x; this implies that, regardless of 
the delay t^, the probability density P(x) is asymmetric. 
By looking at the sign of the noises and the potential 
parameters, one concludes that for r^ = 0a positive drive 
£1 (i.e., tilting the bistable potential V(x) to the right) 
corresponds to rising the right-to-left barrier by a fraction 
proportional to £ 2 = £1, and vice versa for £1 < 0. As 
a consequence the particle is more likely to get trapped 
in the left well and the negative peak oiPJx) overshoots 
the positive one (gating mechanism \l2\ Il9||). 

Two conflicting effects determine the optimal P(x) 
asymmetry: On one side, since the particle takes a finite 
Kramers' time to reach the top of the barrier, delaying £2 
with respect to £1 synchronizes their activation effort thus 
improving the success rate of the noise assisted hopping; 
on the other side, increasing Td larger than r c makes the 
superposition of the two noises statistically incoherent 
and restores the spatial symmetry of P(x). In conclu- 
sion, for low noise intensities, the density asymmetry is 
maximal around Td ~ t c |19|| . 

In the following we investigate the Brownian motion 
on a sinusoidal substrate 

x = a[l + &(*)] cos a; + £i(t) + e(t), (111.2) 

where £1 and £2 are colored noises with correlation time 
t c and relative delay r^, and e(t) is a Gaussian station- 
ary noise with (e(t)) = 0, (e(t)rj(0)) = 0, and (e(i)e(O)) = 
2DS(t). Throughout this section we switch off the ther- 
mal noise e, that is we set D = 0. 

By the same line of reasoning we outlined for the con- 
fined process (IIILlf) . we expect that the Brownian dy- 
namics (fill . 2 1) gets rectified to the left with negative net 
velocity (x(Td)}. Numerical simulation confirms our ex- 
pectations, i.e., (x(Td)} < 0. Figure 0] illustrates the 
dependence of v(rd) = \(x(Td))\ on Td for high [panel 
(a)] and low noise intensities [panel (b)]; in both pan- 
els t c = 0. The curves v(Td) are singular at the origin: 
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FIG. 4: (Color online) Characteristics curve v(Td) = j(i(T"d))| 
for the process l|III.2ft with a = 1, D = 0, t c = 0, and different 
noise intensities (Qi, Q2): (a) high noise levels, resonant tails; 
(b) low noise level, exponential tails versus fitting law lIHI.i 
(dashed lines). In both panels (x(jd)) is negative. 



£2, which occurs only because of the nonlinearity of the 
substrate V(x). This and related aspects of the present 
rectification mechanism will be discussed in Sec. IV. 

On computing the characteristics (or response) func- 
tion v(Td), we remark that the net current (x(rd)) obeys 
obvious symmetry relations that allow us to restrict the 
parameter range to explore: (i) (£(— Td)) = (x(rd)); (ii) 
(x(Td)) — > —(x(Td)) upon changing the relative sign of £1 
and £2- Symmetries (i) and (ii) are easy to prove both 
analytically and numerically (not shown). 

The peculiar Td dependence of the rectification func- 
tion v(rd) is important in view of practical applications. 
Indeed, in many circumstances, it would be extremely 
difficult to recycle a control signal &(t) so that < t c ; 
stated otherwise, measuring v(0) requires a certain de- 
gree of experimental sophistication. On the contrary, if 
we agree to work on the resonant tail of its response curve 
v(rd), a rectification device described by the SDE (|III.2|) 
can be operated with less effort; the net ouput current is 
not the highest, as v(rd > 0) < v(0), but is still apprecia- 
ble and, more importantly, stable against the accidental 
floating of the control parameter Td- 

In this sense the scheme represented in Eq. 



pIL2|l is 

a simple-minded attempt at implementing the operation 
of a Maxwell's daemon: The ideal device we set up is 
intended to gauge the primary random signal £1 at the 
sampling time t and, depending on the sign of each read- 
ing, to lower or raise the gate barriers accordingly at a 
later time t + Td, i.e., open or close the trap door. The 
rectifying power of such a daemon is far from optimal; 
lacking the dexterity of Maxwell's " gate-keeper" [2^ , it 
works "in average" like a "dumb" automaton. 



(i) The rectification speed v(0) is fairly large, in close 
agreement with the predictions of the Fokker-Planck for- 
malism [lj (not shown); (ii) For Td — > 0+ the net speed 
vljd) drops to a much lower non-zero value v(0+) = 
lim Tti _>o K^C^d))!- Moreover, for Td > the curves vfa) 
exhibit a persistent tail that cannot be explained as a 
color effect (we recall that here t c = 0). 

In the low noise regime, Q\ <C 2a and Q2 <§; 1, v(rd) 
decays exponentially with fitting law 



v(Td) = v(0+)e~ 



(III.3) 



whereas at higher noise levels the curves v{jd) develop a 
resonant behavior. The low-noise Td dependence shows 
an obvious resemblance with the Brownian relaxation in a 
parabolic potential driven by a linear superposition of £1 
and £2, see in particular Eq. (|II.3(I : indeed, when sitting 
around the bottom of the substrate wells, the particle 
is subjected to two effective additive noises (harmonic or 
Gaussian approximation), like in Sec. II. In such approx- 
imation the exponential decay (fill . 3|) can be explained 
with the progressive decorrelation of the recycled versus 
the primary noise with increasing Td- The resonant Td 
dependence of the rectification effect at higher Qi,Q2, 
instead, can be traced back to an optimal synchroniza- 
tion of the additive, £1, and multiplicative fluctuations, 




FIG. 5: Characteristics curve v(Td) for the process IIII.2I 
with a = 1, D = 0, Qi = Q2 = 1, and different autocorreal- 
tion times r c . 



IV. RECTIFICATION PROPERTIES 

We discuss now the dependence of the response func- 
tion v(rd) on the remaining parameters of the process 
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FIG. 6: (Color online) (a) Characteristics curve v(Td) for 
the process 1111.21 with a = 1, Qi = Qi = 1, r c = 0, and 
different D. (b) v vs. Q2 with Qi = 1 (empty symbols) 
and vs. Q\ with Q2 — 1 (filled symbols) for different values 
of (jd,T c ): squares (0.2,0.2), circles (0.0,0.2), and triangles 
(0.2,2.0). The remaining simulation parameters are D = 
and a = 1. 

(1) Dependence on the noise autocorrelation time r c . 
As reminded in Ref. [l9j . a finite correlation time t c of 
the cooperating noises £1 and £2 enhances the asymme- 
try of the system response against the relative delay Td- 
In particular, the discontinuity v(0) versus v(0+) in Fig. 
EJa) for r c = 0, is replaced by a continuous drop of v(i~d), 
which approaches its resonant tail down from v(0) on a 
scale Td ~ r c ; conversely, the tails of vfa) decay exponen- 
tially more slowly on increasing r c , until their resonance 
peak disappears completely; 

(2) Dependence on the noise crosscorrelation. So far, 
the recycled and primary noise have been assumed to be 
fully cross-correlated, that is £1 and £2 are both propor- 
tional to the same source 77; the relative time delay Td de- 
termines an effective decorrelation in their action on the 
diffusing particle. We now switch on the thermal noise 
e(t) in Eq 



111.2(1 : The recycled noise £2 turns out to 
be only partially correlated with the total additive noise 



£1 + e, which drives the process with intensity Q\ + D. 
Thermal noise tends to foil the synchronized effort of 
£1 and £2; this is a circumstance of current interest in 
real experiments [T^| . On decreasing the multiplicative- 
additive noise crosscorrelation, e.g., by increasing D, the 
tails of v(rd) in Fig. become less and less persistent, 
i.e., decay faster and faster; their resonance peak shifts 
towards lower Td values, until it merges into the narrow 
peak centered at Td = 23] ; 

(3) Dependence on the recycled-to-primary noise inten- 
sity ratio. The net rectification current speed v exhibits 
a clearcut resonant dependence on both Q\ and Q2, sep- 
arately, regardless of the time constants r c and Td - see 
Fig. Elb). This is certainly true for Td = 0, as known 
from the Fokker-Planck formalism |2fJ, but such an ef- 
fect is more pronounced for intermediate T c ,Td, namely 
< t c , Td < a. Under optimal rectification the maxima 
of v vs. Q in Fig. EJb) occur, as expected, for Qi < 1 
and Q 2 2a. 



V. CONCLUSIONS 



The model <|III.2|) can be regarded as the prototype of a 
rectifying device; inspired by biological systems |2f|. this 
scheme can be employed to design and operate artificial 
rectifiers, for instance, of magnetic vortices and colloidal 
particles. The dependence of the response function v(T d ) 
on the noise parameters, items (l)-(3) of this section, 
indicated how to optimize the net rectification current 
across the device. 

The present investigation was based on the Langevin 
equation formalism - see the analytical predictions of 
Sec. II and the numerical results of Sees. II-IV. Al- 
ternately, one could try to formulate the description of a 
stochastic process x(t) driven by delayed noises in term 
of the probability density P(x,t), namely by generaliz- 
ing the Fokker-Planck (FP) equation formalism [fj. For 
instance, under certain restrictions, the linear SDE ijll.I I) 
corresponds to the time-delayed FP equation 



d x P(x, t) = ad x xP(x, t) + (Qi + < 
±^QiQ 2 d 2 x [P(x, t + T d ) + P{x, t 



h)d 2 x P{x.t) 
-r d )\. (V.f) 



On assuming the asymptotic time-dependence P(x, t) — 
e~ at P(x), for t — ► 00 the FP equation l|V.I|) approaches 
the standard form 



d x P(x, t) = ad x xP(x, t) 
+[(Qi + Q2) ± VQ^h(e aTd 



(V.2) 



]%P(x,t), 



which correctly reproduces the long time behavior of Eq. 
(|II.4|I - but not the stationary quantity Eq. i|II.3|) . A 
rigorous FP equation formalism for SDE with delayed 
noises, valid at all times, is the subject of ongoing re- 
search. 
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